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Sound  generation  in  the  convection  of  turbulence  into  a  flow 
intake  duct  is  discussed  by  examination  of  an  idealised  problem. 
Ideal  fluid  is  in  motion  with  uniform  low,  subsonic  velocity 
above  and  parallel  to  a  plane,  rigid  wall.  A  thin,  rigid, 
semi-infinite  plate  is  parallel  to  the  wall  and  the  fluid  flows 
past  its  leading  edge.  The  radiated  sound  is  calculated  for  a 
turbulent  eddy  modelled  by  a  weak  line  vortex  which  is  allowed  to 
convect  passively  past  the  leading  edge  of  the  plate.  Account  is 
taken  of  the  contribution  to  the  radiated  sound  from  the 
disturbance  produced  by  the  convecting  vortex  in  the  boundary 
layers  on  each  side  of  the  plate  by  means  of  Howe's  (1981)  theory 
of  displacement  thickness  fluctuations,  the  strength  of  the 
disturbance  being  fixed  by  a  leading  edge  Kutta  condition.  It  is 
concluded  that  the  predicted  level  of  the  radiated  sound  is 
substantially  reduced  due  to  the  boundary  layer  disturbances 
relative  to  when  they  are  neglected.  The  case  of  the  convection 
of  a  frozen  two-dimensional  gust  is  also  considered.  Examination 
of  the  analogous  problem  of  plane  wave  radiation  from  the  duct 
shows  that  the  presence  of  displacement  waves  enhanced  both  the 
far-field  intensity  in  the  ambient  fluid  and  the  reflected  field 
within  the  duct. 
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§1  Introduction. 


Sound  produced  when  a  turbulent  eddy  is  ingested  by  a  duct  inflow 
such  as  the  intake  duct  of  an  engine  or  the  entry  of  an  exhaust 
nozzle,  is  of  interest  in  many  industrial  systeas.  Of  concern  in 
similar  circumstances  are  the  sound  fields  which  are  radiated  and 
reflected  at  the  entry  of  a  flow  intake  duct  due  to  acoustic 
sources  within  the  duct,  such  as  a  fan  or  propeller  say.  The  flows  of 
interest  are  often  of  low  Inch  number  and  the  characteristic  acoustic 
wavelength  is  large  relative  to  the  duct  width. 

The  direct  sound  radiation  from  a  compact  turbulent  eddy  convected 
past  a  rigid  body  can  be  thought  of  as  the  sound  which  is  produced  in 
the  absence  of  any  interaction  between  the  eddy  and  the  flow  in  which 
it  is  convected.  It  is  governed  by  the  strain-field  produced  by  the 
presence  of  the  body,  the  details  depending  on  the  rate  of  working  of 
the  eddy  turbulent  stresses  in  that  field  (Howe  1975).  For  Instance 
in  the  two-dimensional  problem  examined  by  Howe  (1975)  in  which  a  line 
vortex  moves  around  the  edge  of  a  thin  semi-infinite  rigid  plate,  no 
sound  is  produced  if  the  vortex  follows  a  potential  flow  streamline 
around  the  edge;  sound  is  only  produced  when  the  vortex  cuts  across 
such  streamlines. 

In  the  presence  of  flow  there  are  additional  sound  sources  which  are 
closely  related  to  the  Incident  vortical  (or  acoustic)  disturbance. 
These  arise  because  of  f low/disturbance  coupling  at  any  edges  of  the 
object  where  vortlclty  can  be  generated  and  ejected  into  the  flow,  a 
process  which  is  essentially  viscous  controlled.  The  mechanism 
provides  a  means  of  transferring  energy  that  is  well  illustrated  in 
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the  impingement  of  a  sound  wave  -on  the  leading/ trailing  edge  of  a  thin 
seal-infinite  rigid  plate  in  a  low  subsonic,  grazing  aean  flow  and 
respectively  results  In  a  net  production/absorptlon  of  acoustic 
energy,  Howe  <1981 ). 

In  this  paper  we  discuss  sound  production  by  the  Interaction  of  a 
vortical/acoustic  disturbance  with  an  otherwise  uni fora  duct  inflow. 
Though  substantial  aaounts  of  acoustic  energy  can  be  produced  when  the 
incident  flow  is  not  self-preserving,  as  in  the  jet-edge  tone 
configuration  (Crlghton  1984),  in  the  cases  examined  here  the 
disturbances  are  assumed  to  be  sufficiently  weak  that  both  the 
disturbance  and  acoustic  particle  velocity  are  small  relative  to  that 
of  the  mean  flow. 

The  production  of  sound  in  the  passive  convection  of  a  line  vortex 
by  a  low  Inch  number  mean  flow,  of  speed  U,  past  the  mouth  of  a  rigid, 
two-dimensional,  semi-infinite  duct  is  examined  in  82.  Longitudinal 
standing  waves  are  avoided  in  this  first  analysis  by  considering  a 
duct  of  infinite  length.  A  more  realistic  model  will  also  require  a 
three-dimensional  treatment  to  obtain  predictions  of  the  sound 
produced  in  the  ingestion  of  a  general  turbulent  field  by  a  duct 
Inflow.  The  duct  is  formed  by  an  infinite  plane  rigid  wall  and  a  thin 
semi-infinite  rigid  plate,  parallel  to  the  wall  and  a  distance  H  above 
it  in  the  fluid.  H  is  assumed  small  relative  to  any  relevant  acoustic 
wavelength.  The  mean  flow  lach  number  I  =  D/c  <c  is  the  sound  speed) 
Is  small  and, for  convenience  in  the  analysis; terms  of  order  P  are 
neglected  relative  to  unity.  Thus  c  and  the  aean  fluid  density  p  can 
be  taken  as  constant.  The  theory  presented  is  linear  with  respect  to 
the  acoustic  particle  velocity. 
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The  flow  model  Is  introduced  is  S2.1  and  the  analytical  problea  Is 
formulated  in  S2.2.  In  S2. 3  Howe's  <1981)  theory  of  displace sent 
thickness  fluctuations  is  applied  to  deter  sine  the  effect  on  the 
radiated  sound  field  of  disturbances  propagating  in  the  plate  boundary 
layers.  Physically  such  disturbances  arise  fros  vorticity  production 
at  the  plate  leading  edge  by  the  action  of  viscosity  in  response  to 
the  unsteady  velocity  field  of  the  passing  vortex.  The  wall  boundary 
layer  is  presumably  well-established  and  stable,  unlike  the  thin 
developing  plate  bounds  -f  ayers  which  are  linearly  unstable.  In  the 
real  flow,  acoustic  and  the  ml  boundary  layers,  analogous  to  a 
Stoke' g  layer,  will  also  be  present  on  the  plate  and  wall  (due  to  the 
need  to  satisfy  the  no-slip  condition  at  their  respective  surfaces! 
and  will  be  extresely  thin  relative  to  an  acoustic  wavelength.  Any 
Interactions  between  vortex  or  sound  shear  waves,  say,  and  boundary 
layer  Reynold's  stresses  will  be  insignificant  over  the  length  and 
tlse  scales  which  characterise  the  Interaction  of  the  vortex  with  the 
duct  south  and  are  neglected  here.  Similarly  aerodynamic  sound 
generation  (Llghthlll  1952)  by  SBan  or  perturbed  Reynold's  stresses  is 
not  Included.  Thus  it  is  not  necessary  to  account  for  the  boundary 
layer  on  the  wall  and  the  predominant  effect  of  the  plate  boundary 
layer  is  accounted  for  by  Howe's  (1981)  displacement  wave  theory. 

In  Howe's  theory,  fluctuations  in  the  displacement  thickness  of  the 
boundary  layers  model  the  affect  of  the  disturbances  in  the  exterior 
potential  flow.  Here,  as  la  Howe  (1981),  the  boundary  layers  are 
modelled  as  uniform  wall  flows  of  speed  T  <  H  and  the  fluctuations  are 
modelled  as  neutrally  stable  novas  of  constant  fora  and  of  long 
wavelength  relative  to  the  characteristic  boundary  layer  thickness. 

The  displacement  waves  carry  perturbation  energy  downstream  into  flow 
regions  where  they  are  dissipated  by  heat,  viscosity  or  more 


realistically  perhaps,  breakdown  into  saall-acale  turbulence.  Thus 
the  effective  Interaction  region  Bust  be  close  to  the  leading  edge  of 
the  plate.  This  physical  assumption  is  incorporated  into  the  analysis 
by  allowing  the  wave-number  of  the  displace  sent  waves  a  snail 
imaginary  part  where  necessary  to  ensure  the  convergence  of  certain 
Integrals.  The  strength  of  the  displacement  waves,  proportional  to 
that  of  the  vortex,  is  fixed  by  the  Kutta  condition  that  the  fluid 
velocity  should  resain  finite  at  the  plate  leading  edge.  Viscous 
effects  enter  the  analysis  only  in  this  indirect  fashion. 

Though  the  leading  edge  Kutta  condition  in  unsteady  flows  has  been 

exaained  by  Goldstein  (1981,1983)  the  issue  of  its  validity  has  not  been 

resolved.  However  Howe's  (1981)  application  of  the  displacement  wave 

theory  incorporating  the  leading  edge  Kutta  condition  to  the  Jet-drive 

aechanisa  of  the  flue  organ  pipe  gives  encouraging  agree went  with 

experiment.  The  edges  of  a  real  duct  mouth  will  have  finite  thickness 

relative  to  the  displacement  wave  length  scale  and  may  be  rounded 

rather  than  sharp,  levertheless  the  mechanism  controlling  the 

flow/acoustic  interaction  at  the  mouth  will  be  essentially  the  same  as 

that  involving  a  thin  plate.  Goldstein  (1984)  shows  that  when  a 

laminar  flow  separates  from  a  smooth  body  surface  and  is  subject  to  an 

external  unsteady  'forcing'  (e.g.,from  a  harmonic  source),  instability 

wives  which  propagate  downstream  are  generated  and  are  coupled  to  the 

source  at  tha  separation  point  by  viscous  effects.  In  particular  for 

a  harmonic  source  the  Kutta  condition  is  satisfied  at  the  separation 

1/4 

point  provided  that  1  <<  S  <(8e  ,  where  S,  Ke  are  respectively  the 
Strouhal  and  Reynold's  numbers  based  on  the  streamwlse  body  length 
scale  and  the  upstream  flow  velocity,  i  comprehensive  review  of  the 
Kutta  condition  la  unsteady  flows  is  provided  by  Crlghton  (1985). 
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The  solution  for  a  fourier  tlse  transform  of  the  sound  field  Is 
obtained  in  S2.4  by  use  of  a  low  frequency  Green's  function  and  is 
found  to  have  two  components;  the  direct  radiation  froa  the  vortex  and 
the  contribution  froa  the  displacement  waves  an  the  plate  boundary 
layers.  On  application  of  the  futta  condition  we  find  that  the  total 
field  is  siaply  the  direct  radiation  multiplied  by  a  factor  equal  to 
(l-V/U*)  =  .4,  for  a  characteristic  value  of  V/De  where  0.  (=  D>  is 
the  convection  speed  of  the  vortex.  Thus  the  presence  of  displacement 
waves  substantially  reduces  the  radiated  field.  Comparison  of  the 
direct  radiation  is  made  with  the  solution  obtained  by  Cannell  and 
Ffowcs-Villiams  <1973)  in  the  absence  of  mean  flow  for  a  line  vortex 
exhausting  out  of  the  duct.  Explicit  forms  of  the  far-field  due  to 
the  line  vortex,  both  in  the  duct  and  outside  it,  are  obtained  in  52. 5 
and  are  valid  when  the  duct  is  wide  relative  to  the  displacement  wave 
length  scale,  characterised  by  the  minimum  distance  of  the  vortex  from 
the  plate.  This  shows  that  it  makes  little  difference  whether  the 
vortex  moves  above  the  plate  or  beneath  it  into  the  duct. 

The  results  of  £2.5  are  used  in  £2.6  where  the  radiated  field 
produced  in  the  convection  into  the  duct  of  an  incompressible,  frozen 
two-dimensional  gust  is  discussed.  The  flux  of  acoustic  energy,  both 
to  the  far-field  in  the  duct  and  outside  it,  is  calculated  in  £2.6. 

For  the  line  vortex  case  these  fluxes  are  compared;  the  former  is 
proportional  to  I,  the  latter  to  IP. 

In  £3  a  low  frequency  plane  sound  wave  propagating  out  of  the  duct 
replaces  the  line  vortex  as  the  'incident'  sou  id  source.  The 
displacement  wave  theory  outlined  above  is  employed,  la  this  case  the 
presence  of  displacement  waves  results  in  an  increase  in  both  the  far- 
field  intensity  outside  the  duct  and  the  amplitude  of  the  reflected 
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field  within  the  duct.  There  is  experimental  evidence  (Davies  1987) 
of  such  an  increase  in  the  reflection  coefficient  in  plane  wave 
reflection  at  a  flow  Intake.  The  calculation  of  S3. 2  reveals  that  the 
net  acoustic  energy  flux  away  iron  the  plate  leading  edge  is  positive; 
acoustic  energy  is  generated  by  the  displacement  thickness  waves. 

This  Is  in  accord  with  the  results  of  Howe  (1981)  for  the  interaction 
of  a  plane  wave  with  a  plate  leading  edge  and  of  Quinn  A  Howe  (1986). 
Appendices  A1-A4  contain  various  analytical  results. 


S2.1  The  flow  model 


The  problem  to  be  considered  is  illustrated  schematically  in 
figure  1.  A  uniform  mean  flow  of  low,  subsonic  velocity  D  proceeds 
along  a  rigid  wall  (at  *2  =  0  in  the  (11,12)  co-ordinate  system  of  the 
figure)  and  encounters  a  two-dimensional  duct  formed  by  a  thin,  rigid 
semi-infinite  plate  at  12  -  B  above  the  wall.  A  turbulent  eddy 
modelled  by  a  line  vortex  of  strength  y/2x  convects  at  speed  0e  along 
a  path  12  =  h.  Since  the  flow  Kach  number  satisfies  K  =  D/c  <<1,  the 
speed  of  sound  in  the  fluid,  c,  may  be  assumed  constant  and  in  tke 
subsequent  analysis  terms  of  order  IP  are  neglected  relative  to  unity. 

The  aim  of  the  analysis  is  to  provide  estimates  of  the  sound 
produced  by  the  interaction  of  the  vortex  with  the  duct  mouth. 

Besides  the  'direct  radiation1  from  the  vortex  there  are  acoustic 
sources  closely  related  to  the  vortex  which  arise  when,  in  response  to 
the  velocity  field  of  the  passing  vortex,  additional  vortlclty  is 
ejected  into  the  plate  boundary  layers  at  the  plate  leading  edge.  The 
assumed  region  of  this  interaction  is  indicated  schematically  by  the 
wavy  lines  in  figure  1.  The  flow  region  of  importance  to  the  acoustic 
far-fleld  is  assumed  to  be  the  vicinity  of  the  duct  mouth,  lx>  <  B, 
and  should  be  confirmed  by  the  results  of  the  analysis. 


The  turbulent  boundary  layers  on  the  wall  and  plate  will  be 
extremely  thin  relative  to  the  acoustic  wavelengths  concerned  and 
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Figure  1.  Schematic  Illustration  of  convection  of  a  line  vortex  by  a 
uniform  mean  flow  into  a  semi-infinite  duct. 
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Figure  2.  Sketch  of  the  idealized  boundary  layer  model  with  a 
displacement  wave. 
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their  interactions  with  the  acoustic  fields  produced  are  assuned  to  be 
incoherent  in  respect  to  the  acoustic  tlae  and  length  scales  of 
Interest.  The  line  vortex  is  assumed  to  be  sufficiently  weak,  compact 
and  Incoherent  with  regard  to  the  turbulent  boundary  layers  that 
additional  sound  fields  arising  froa  perturbed  Beynolds  stresses  can 
be  neglected.  Thus  the  wall  boundary  layer  is  excluded  and  the  plate 
boundary  layers  enter  the  model  only  in  relation  to  the  leading  edge 
interaction  outlined  above. 

The  fluid  Mechanics  governing  the  vortex  is  of  interest  only  in 
the  effective  interaction  region,  1*1  <  H  where  the  vortex  (or  gust  of 
S2.5)  is  assuned  to  convect  passively  at  a  speed  0« .  Variations  in  Uc 
and  curvture  in  the  vortex  path  during  the  passage  of  the  vortex  past 
the  duct  nouth  due  to  an  inage  system  in  the  plate  will  be  snail 
provided  that  y/IH-hl  <<  and  are  neglected.  A  physically  realistic 
convected  disturbance  might  have  a  convection  velocity  of  .50  to  .8U. 
The  simplistic  modelling  of  the  vortex  means  that  the  energetics  of 
the  Interaction,  being  of  second  order  in  the  acoustic  fluctuations 
produced,  is  not  included  explicitly  in  the  flow  model.  Kadiated 
acoustic  energy  produced  in  the  interaction  is  balanced  by  a  decrease 
in  mean  flow  energy  or  energy  of  the  convected  disturbance  which  in 
the  real  flow  might  be  expected  to  be  apparent  in  higher  turbulence 
levels  in  the  plate  boundary  layers  downstream  of  the  edge. 

In  addition  to  X  <<  1  further  physically  acceptable  but  limiting 
assunptlons  enter  the  analysis.  Since  typically  turbulence  levels 
will  not  exceed  10%  of  the  mean  flow  velocity  we  take  uh  <<  D,  uh 
being  the  characteristic  hydrodynamic  disturbance  velocity.  The 
acoustic  particle  velocity,  IV# I  <#  is  a  perturbation  potential)  will 
also  be  small  relative  to  0  and  0*  =  0(0).  These  Inequalities  may  be 
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expressed  as  t^fl,  uh  <<  0  <<  c.  The  neglect  of  teras  of  order  K2 
relative  to  unity  Is  for  mathematical  convenience.  The  assumption  of 
an  acoustically  compact  Interaction  region  should  be  satisfied 

I 

provided  that  the  duct  height  Is  saall  relative  to  any  relevant  { 

acoustic  wavelength,  i.e.  H  <<  c/o  for  a  characteristic  radian 

frequency  a.  Effects  due  to  theraal  diffusion  processes  can  be 

assuaed  to  be  insignificant  in  view  of  the  relatively  short  tiae 

envisaged  for  the  passage  of  the  disturbance  through  the  region 

I  *1  <H. 


§2.2  The  analytical  problem 

Given  this  modelling  we  consider  the  analytical  problem  illustrated 
schematically  in  figure  lj  ideal  fluid  is  in  motion  in  the  positive  x, 
direction  of  the  (Xi.Xa)  co-ordinate  system  with  unifora  low,  subsonic 
velocity  U  above  a  rigid  wall  at  xz  -  0.  A  thin  rigid  semi-infinite 
plate  is  at  X2  =  H,  Xi  >0  and  is  parallel  to  the  wall.  A  line  vortex 
of  strength  t/2x  convects  at  a  constant  speed  Ut  at  a  constant 
distance  h  above  the  wall. 

Vhen  visco-thermal  effects  are  neglected  the  stagnation  enthalpy 
B<x,t),  defined  by 


B  =  w  +  v2/2  (2. 1) 

where  w  is  the  specific  enthalpy  and  2.  is  the  fluid  velocity,, will 
satisfy  the  Inhomogeneous  convected  wave  equation  (Howe, 1975); 
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B<x,  t '  =  -div(u  x  e). 


(2.2) 


'V  +  V  - 

(  1^  +  1&  >2 

.  ^Xi2  ^Xa2 

t  cit  ixi)  j 

B  Is  exploded  here  as  the  fundamental  acoustic  variable.  In 
equation  (2.2)  u  =  curia  is  the  vorticity,  c  Is  the  speed  of  sound 
and  any  be  assumed  to  be  constant  since  the  Xach  number  X  =  U/c  of 
the  flow  satisfies  X2  <<1.  In  irrotational  flow  regions  we  may  take 

B  =  -  ^*At  (2.3) 

where  i  is  a  perturbation  potential  of  the  flow. 

The  source  term  in  (2.2)  will  be  non-zero  only  in  vortical  flow 
regions,  i.e.  at  the  vortex  core  given  by  it  =  Uc(t-to),  Xa  -  h 
(where  to  is  a  constant)  and  possibly  at  the  surfaces  of  the  plate 
x2  =  H±e,  where  additional  vorticity  may  be  introduced  into  the  flow 
at  the  leading  edge  by  the  passing  line  vortex.  The  vorticity  of  the 
convecting  line  vortex,  ui,  say  is 

Hi  -  I  f(  Xi  +  Uc to  -  Oct)  S(X2-h)  a  (2.4a) 


where  a  is  a  unit  vector  out  of  the  plane  of  the  paper  in  figure  1  and 
S  is  the  Dirac  delta  function.  Ui  can  be  expressed  as  a  superposition 
of  vorticity  waves: 


Ut (X.t) 


■f 


Ui  (x,  t,u)du 


where 


(2.4b) 


Ui  =  (Y/2xUc)<(X2-h>exp(iutto-ii+Xi/0c))a  . 


Henceforth  hatted  quantities  will  denote  variables  in  the  frequency 
dam  in  and  the  tlm  dependence  exp(-iut)  will  not  be  denoted 
explicitly.  The  solution  for  the  stagnation  enthalpy,  B  can  be 
obtained  by  superposition  of  the  solution,  B  say,  associated  with  the 
vortical  field  Ut  of  (2.4b).  Ve  note  that 

«i  x  x  =  (y/2x)4 (X2-h)exp(iot to+xi /Del  )n  ,  (2.5) 

where  n  is  a  unit  vector  in  the  positive  xa  direction  of  figure  1. 

S2.3  Boundary  Conditions 

The  far-field  radiated  sound  will  be  calculated  from  equation  (2.2) 
with  the  incident  vorticity  field  Ui  of  (2.4b)  and  appropriate 
boundary  conditions  at  the  plate  and  wall.  On  the  rigid  wall,  x2  =  0, 
the  zero  normal  velocity  condition  requires  that  'bli/'bxa  =  0,  since 

/V  A 

from  (2.3)  and  (2.4b),  B  =  iul.  The  same  condition  applied  on  the 
plate,  Xi  >  0,  X2  =  H  would  lead  to  singularities  in  the  velocity  and 
pressure  at  the  plate  leading  edge,  x  =  (0,H).  The  singularities 
result  from  the  idealization  that  viscous  effects  are  negligible 
whereas  In  practice  viscosity  acts  to  inhibit  such  singular  behaviour. 
The  singularities  are  avoided  here  by  use  of  Howe' 6  (1981)  theory  of 
displacement  thickness  waves.  In  that  theory  the  effect  on  the 
exterior  potential  flow  of  large  scale  boundary  layer  disturbances 
(produced  by  vorticity  generated  at  the  leading  edge)  is  modelled  by 
fluctuations  in  the  displacement  thickness  of  the  boundary  layers 
emanating  from  the  edge.  Greater  detail  of  the  theory  is  given  by 
Howe  (1981)  where  an  Idealized  model  of  the  boundary  layer,  sketched 
in  figure  2  and  adopted  here,  is  employed.  In  the  model  a  boundary 


layer  wave  of  the  fora  Aexp(lxxi),  y  =  x (u) ,  propagates  above  a 
uni fora  wall  flow  of  speed  V  with  V  (  0,  0  being  the  speed  of  the 
exterior  uni  fora  flow.  Howe  shows  that  though  a  wide  range  of  values 
of  x(u)  are  possible  for  each  o,  as  )[{  i  0  (6  is  the  characteristic 
boundary  layer  thickness!  the  various  aodes  coalesce  and  x  "*  "  /?•  We 
take 

X  =  u  /  T  (2.7) 

and  suppose  that  V  =  0.60,  a  value  suggested  by  the  work  of  Bull 
(1967)  and  Blake  (1970).  Thus  the  boundary  condition  at  control 
surfaces  defined  by  x2  -  Htd ,  z,  >  0,  Just  outside  the  thin  boundary 
layers  on  the  plate  is 

iB/ix2  =  loi'(o)exp(ixXt ),  Xi  >  0,  Xa  =  H±0,  (2.8) 

where  J/(u) ,  a  weasu re  of  the  strength  of  the  disturbances  in  the 
boundary  layers,  is  to  be  fixed  by  the  Kutta  condition  that  the  fluid 
velocity  is  finite  at  the  leading  edge.  In  (2.8)  it  has  been  assuned 
that  the  condition  can  be  applied  on  the  surface  of  the  plate. 

S  2.4  Tike  solution  for  a  vortex  wave 

Introduce  the  degenerate  Prandtl-Glauert  transforation  (in  which  IP 
is  neglected  relative  to  unity) 

W  A 

B(x>  =  B(*)exp(ioDxi /ca> ,  (2.9) 

The  governing  wave  equation  (2.2)  becoaes 


(2. 10) 


f  **  +^2  +  f  ' 

1  ixi2  ixa2  c2  , 


ll(x>  =  -div(£  x  z>exp(luDxi /c2) 


where  again  V2  has  been  neglected  relative  to  unity.  Ve  define  G(x;y> 
as  the  solution  of 


+  i2  to2' 
.  iyi2  iya2  c2  . 


G(X;y) 


-6 (x-y) . 


(2.11) 


satisfying 


ic/iy2  =  0,  ya  =  0,  lyi  I  <  •  , 


^G/iy*  =  0,  ya  =  H,  yi  >  0  , 


(2.12) 


and  the  radiation  condition  that  G(x;y)  should  have  outgoing  wave 
behaviour  in  y  space  for  olyl/c  >>  1,  and  that  as  yi  -*  •  with 
lyal  (H,  G(x;y>  -*  constant. exp <luyi/c>.  G(x;y)  is  calculated  in 
Appendix  A2. 


From  equations  (2.10),  (2.11)  we  find 


tf<x>  =j(G<x;y)V2E(y>  -  B(y)V2G(x;y)  +  G(x:y)div(uxy)exp(iuDy,/c2))  dv. 


(2. 13) 


The  integration  in  (2. 13)  is  over  the  region  of  fluid  V'  which  is 
enclosed  by  a  control  surface  S.  S  consists  of  the  rigid  wall  ya  =  0, 
a  surface  y,  -  constant  =  ■  for  I yal  <  H,  control  surfaces  at  ya  =  H±£ 


Figure  3. 


Sketch  of  the  control  surfaces  of  equation  (2.14)  and  the 


vortex  wave. 
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for  j\  >  0  (around  the  plate)  and  a  circular  arc  centred  on  (0,H>  to 
give  a  closed  circuit.  Application  of  Green's  theorea  to  the  fluid 
within  S  gives  froa  (2.13) 


£<X> 


=  Jtt.  <G<x; 


y)VB(y)  -  B(y)7  G(xij)  +  G(x;  Jt>  (uxs)exp(loDy,  /c2)  )dS 
V (G(x;y)exp(iuyi/c2))  dV  ,  (2.14) 


-^<U  x  x) . 1 


where  (see  figure  3)  a  Is  the  unit  noraal  to  S  In  the  direction  of  the 
interior  of  S.  loting  that  vortical  fluid  enclosed  by  S  is  that  of 
the  vortlclty  wave,  we  find  on  using  the  radiation  conditions  on 
G(x;y>  and  the  conditions  (2.12)  that 


B(x)=  -{(G(x;j)liB(y) 
^ya 


dy,  -yexp(iuto)[ \c  (x;  y)exp(iuyi  (l+lIB«/c)  /0c  >dyi  . 


2 x 


y* 


(2. 15) 

where  continuity  in  )>B/iy2  across  the  plate  (c.f.  (2.8))  has  been 
used.  In  this  equation  the  first  integrand  is  evaluated  at  ya  -  B, 
the  second  at  y2  =  h.  Ve  have  used  equation  (2.5)  for  the  tera 
(u  x  v)  and  the  notation 

4 

if)  =  f(y.,H+0)  -  f(y,,H-0). 


lb 


(2. 16) 


The  plate  boundary  condition  <2. 8)  is  substituted  In  (2.15)  and  the 
Prandtl-Clauert  transformation  (2.9)  reversed  to  give  the  actual 
enthalpy  as 


£<x>  = 


Bd(x>  +  Bw(x>  i 


<2. 17  > 


where 


Bw<x>  =  -texp<io<to-0ci /c3)  )/exp<  (iuyi  )t  l+KD./cl /0«  Ag(x;  yi .  h)dy>  , 


2c 


(2.1B> 


i 


B« <x>  =  -ioVe*p(-luDxi/c3)  ecp<lxyi+iuOy>/ca>tG<x;y)l  dyi 


(2. 19) 


Bw<X>  Is  the  direct  contribution  to  the  sound  field  from  the 
vortlclty  wave  of  <2. 4b)  while  iL(x)  is  due  to  the  disturbed  boundary 
layers  on  the  plate.  The  Integrals  of  (2.18),  (2.19)  are  approximated 
by  use  of  equation  (42.9)  of  ippendlx  42  which  Is  valid  for  positions 
X  in  the  far-fleld,  y  In  the  near-fleld  and  gives 

(G<x;y>l  «  o(x,u)l£(j)t  , 


<>G(x;y>/  ya  *  a(c,u>^e<J>/iya  , 


(2. 20a, b> 


I* 


for  ulxl/c  >>  0(1>,  ulyl/c  «  0(1).  "J  Is  a  potential  function  of  an 
Inflow  Into  the  duct  (unifora  when  ly>l<H,  y <  -*•  )  end  is  discussed  in 
Appendix  Al.  a  is  the  function 


with 


10  jo 

a(j, u)  =  -u(Ho(ulx>  /c)  +  l  <x, /I  x>  >Hi  (ul  xl  /c)l  , 


2c  (1+iuHf (u)/cc) 


f(o)  =  .5772-ix/2+ln(uB/2xc). 


(2.21a) 


(2. 21b) 


<«  <« 

(c.f  (A2.8), (A2. 12)  of  Appendix  A2).  Bo,  H<  are  Hankel  functions  of 

(.1 

the  first  kind.  The  tera  Involving  Bo  Is  a  aonopole  scattered  field 
and  that  involving  tfi*  is  a  dipole  scattered  field. 


(1)  To  evaluate  Bo(x>:  the  direct  radiation 

Equations  (2.20b)  and  (2.19)  give  (with  results  (A1.4),  (A1.6)  oi 
Appendix  Al) 


B*<z)  =  (T/2x>y{a(x,u)exp<iu(to-l(x,/c) )  , 


where  (c.f  definition  (A4.9)  of  Appendix  A4) 


Jexp(l*iy,/U.)Ie(l/(l-Zoa>)dy,  , 


(2.22a) 


(2. 22b> 


and  Zo(y<)  is  the  iaage  of  the  point  y«+Jh,  J= (-1 )“,  under  the 
tranef orantlon  (A1.2)  of  Appendix  Al.  %*(X>  corresponds  to  the 
radiated  field  calculeted  by  Canaell  and  Ffowca-Vllliaas  (1973),  by 


■cans  of  Batched  asyaptotic  expansions,  for  the  case  when  a  vortex 
exhausts  froa  the  duct  in  the  absence  of  aean  flow,  its  action  due  to 


i wages  in  the  duct  walls.  Equation  (2.22a)  agrees  with  their  results 
when  account  is  taken  of  the  different  vortex  actions,  save  that  the 
tera  f (o)  in  a(x,  u>  of  (2.21)  is  absent  In  their  calculation.  The 
discrepancy  appears  to  arise  because  of  their  erroneous  exclusion  of  a 
pressure  variation  which  is  unlfora  over  the  whole  duct  south  (c.f. 
equation  (3.15)  of  Cannell  and  Ff owes- Vi 11 Ians  (1973)). 


(11)  To  evaluate  B_(x>:  the  dlsplaceaent  wave  contlbutlon 


■eglectlng  IP  relative  to  unity,  we  find  froa  (2,19),  (2.20a)  that 


K_(x> 


f* 

-luv (u)exp(-luKx, /c)a(x.u)  expdxy i  >1 nl_  dy> 


(/ 


lwp(u)exp(-lulxi /c)a(x.u)  I  f  d(exp(lxy<)>  df 

J 


(2.23) 


where  now  the  integration  variable  is  the  potential  f.  The 
dlsplaceaent  waves  Bust  decay  as  y,  -*  •  due  to  dissipation  by  heat  and 
vlacoalty,  or  breakdown  Into  turbulence,  and  this  is  aodelled  here  by 
allowing  u  a  aaall  positive  iaaglaary  part  (for  u>0>  in  the  evaluation 
of  the  integral  of  (2.23)  to  give 


*.(X> 


-<wV(w>/x>exp(-iufe,  /c)a(x.u)  fexpdyy,  )d£ 


f 

<x,«>  • 


Caspar i son  with  (14. 5)  of  Appendix  14  leads  to  the  result 

B»<X)  =  -{uV(o)H/x*)exp(-luhxi /c)a(x,u)exp[L(-ixH/*))  .  (2.24) 

where 

L(x>  =  ln(T(x)l  -xlnx  +x  ,  (2.25) 

and  T  is  the  gasan  function  (Cradshteyn  and  Kyxhik  1980, p933) . 

The  calculation  of  Ippendlx  14  fixes  V (u)  by  application  of  the 
Kutta  condition  that  the  fluid  velocity  resale  finite  at  the  plate 
leading  edge  and  gives  (c. f. equation  (14.10)) 

V<o)  =  «(l/2*)expl  iuto-L('ixH/*)]){/HUc  .  (2.26) 

On  substitution  in  (2.24)  we  find 

B~  <X>  =  -  (»/0«  )  (t/2x))/a(x  ,  u)exp(iu(  to-Kxi  /c)  ,  (2.27) 

which  is  valid  in  the  far-fleld,  ul*l/c  )>1.  Cosparlson  with  (2.22) 

reveals  that 

B(x>  =  Bu<x)  t  Bw(x)  =  (1-y/Uc  )Bu(x)  .  (2.28) 

The  ters  coataiaiag  V  arises  froa  the  displacesent  waves  so  their 
presence  has  the  effect  of  reducing  the  level  of  the  radiated  sound  by 
a  factor  of  (1-1/0.)  which  for  a  characteristic  value  of  I  -  0.60  and 
0«  *  0  la  A '10.  The  reduction  ia  greater  if  0.  <  0  and  for  T  =  Dc  no 
aouad  is  radiated,  as  if  the  plate  were  absent.  It  is  interesting  to 
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compare  this  result  with  those  of  Howe  (1976)  for  the  case  when  a  line 
vortex  convects  past  the  trailing  edge  of  a  thin  seal-infinite  rigid 
plate  in  a  low  subsonic  grazing  aean  flow.  Application  of  the 
trailing  edge  Kutta  condition  required  vorticity  to  be  shed  from  the 
edge  and  to  convect  at  the  aean  flow  velocity  in  a  vortex  sheet  wake 
downstreaa  froa  the  edge.  Howe  found  that  no  sound  was  radiated  if 
the  convection  velocity  of  the  vortex  and  shed  vorticlty  were  equal. 


S2.5  The  solution  for  a  convected  line  vortex  upstreaa  disturbance 


The  stagnation  enthalpy  B(x, t)  for  the  line  vortex  of  (2.«a>  is 
obtained  by  integration  of  B<x>  over  «.  Ve  suppose  that  the  duct  is 
extreaely  wide  on  the  hydrodynamic  length  scale  D/o,  in  which  case  the 
integral  of  (2.26b)  is  approximated  by  equation  <A4.14>  of 
Appendix  A4.  Ve  find  froa  (2.22),  (2.20)  that 


B(x,  t)s  (1-7/0.  )  (y/4x)  <H0«)“(l-i>/  a(x,u)exp(-iwt>(  exp(-oHX/Bt )  J  Ido  tc.c 


lul** 


(2.29) 


j 


where 


t  =  t-totlXi/C, 


and 


X.  =  ll-h/HI,  X-  =  (1+h/H) . 


(2.30) 


The  difference  in  the  Integrand  is  to  be  taken  between  X=I*,  c.c 
denotes  the  complex  conjugate  of  the  preceding  expression  and  we  have 


i 


used  «*(x,u)=a(x,-u>.  K,*(o)=|£(-u)  <c.f.<2.21)  and  (2.22b>>.  An 
asterisk  denotes  the  complex  conjugate. 

(1)  Far-fleld  outside  the  duct 


On  substitution  for  ct<x, u)  from  (2.21)  into  (2.29)  we  obtain 


B(x,t) 


(1-V/0*  ItEc-d+cose) 


4x (B/Uc ) (xl xl /H)“ 


1  exp (-1st) (exp(-sX» )-exp(-sX->) Ids+c. c. 

J  —  (1+ilUf  (sOc/H)/x) 


(2.31) 


The  asymptotic  form  of  the  Hankel  functions  has  been  used  (Ibraaowi tz 
and  Stegun  1964, p364).  In  (2.31)  we  have  introduced  IL ,  8  and  ?; 

I U  =  0*/c,  8  =  arcos(x>/lxl  )  ,  T  =  Ddtl/H,  (2.32a) 

where  since  K3  <<  1 

It]  =  t-to-lxt/c  +Kxi/c  =  t-to-lxl/c(l+Icos(8)>.  (2.32b) 

f  Is  defined  by  (2.21b).  Owing  to  the  exponential  decay,  the  major 
contribution  to  the  integral  in  (2.31)  is  from  s  <<  1/1.,  where  Is 
the  minimum  distance  of  the  vortex  from  the  plate  as  a  fraction  of  the 
duct  height.  Correct  to  neglect  of  terms  of  order  (Ks)3la(Bs)  we  find 
after  some  calculation  and  use  of  Gradshteyn  and  Kyzhlk  (1980,  p573) 
that 

B(x,t>  =  le(F(T,I*)-F(T,I-))  . 

Aw  =  (l-?/Ue>  C|/2x)I«*(l+coe8) . 

(xlxl/H)*(H/0.) 

XX 


(2.33) 


Here 


F(t,x)  =  t  -lUd  fxg(t,x)-tln<X./2x(x2+t;*)“>1  , 


x2+t2  x  dti  xa+ta 


and 

g (x,  t )  =  x/2  -  arcos<x/  (x2+t2)**) . 


<2. 34a) 


<2. 34b) 


Since  T  and  X*  are  non-dimensional  variables,  so  is  F  In  (2.33)  and 
An  Is  of  order  u0,  where  u  is  the  characteristic  particle  velocity  of 
the  vortex.  The  paraaeter  Y  =  De(t]/H  is  the  retarded  streamwise 
position  of  the  vortex  as  a  fraction  of  the  duct  height.  The  tern 
F(Y,X»)  arises  from  diffraction  by  the  plate  and  dominates  the  term 
F(Y,X->  (due  to  the  presence  of  the  wall  at  xs  =  0)  when  the  vortex 
passes  close  to  the  plate.  B/An  describes  the  dependence  of  the  far- 
field  enthalpy  on  the  retarded  vortex  positon.  As  III  ^  ®  and  the 
vortex  is  far  from  the  leading  edge  the  enthalpy  decays  with 
B/An  =  <4h/H) /III3.  Since  (c.f(2.3Aa)) 


F(t.x)  *  {t<l+m/x)-JM  l+ln(2xx2/JU>l/x)/xa  ,  t  <<  x, 


when  D-(tl  <<  IH-hl  and  the  vortex  approaches  its  minimum  distance 
from  the  plate  edge,  B/An  becomes  small  with 

B/An  =  (OOMtl/IH-hl  )+0<X))/l  1-h/HI . 

B/An  will  have  a  maxlnua/minimua  when  the  vortex  is  respectively 
downstream/upstream  of  the  plate  edge  and  for  X  <<  1,  I 1-h/HI  <<  1 
this  occurs  at  Y  =  ±1 1-h/HI  where  B/An  =  ±1/21 1-h/HI.  lote  that  the 
assumption  that  the  vortex  convects  approximately  at  the  mean  flow 
velocity  relies  on  the  assumption  Y/H,  f/IH-hl  <<  0,  so  that  the  self- 
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Induced  action  of  the  vortex  is  negligible.  Consequently  the  results 
presented  here  are  not  valid  as  h-tO  for  fixed  Y  and  0. 

The  characteristics  of  B/A«  discussed  above  are  evident  in 
figure  4(a,b)  where  respectively  He  =  . 1,  .001  and  B /A«  is  plotted  as 
a  function  of  Y  =  0*1  tJ/H,  I7I<1,  for  various  values  of  h/H. 

Figure  4<a)  shows  that  it  Bakes  little  difference  whether  the  vortex 
passes  above  or  below  the  plate  though  the  amplitude  IB/AkI  is 
slightly  greater  in  the  foraer  case.  Figure  4(b)  indicates  that  for 
He  <<  1,  IBI  does  not  depend  on  whether  the  retarded  vortex  position 
is  upstream  or  downstream  of  the  duct  mouth.  B/A«  is  plotted,  again 
as  a  function  of  Y  in  figure  5  for  fixed  h/H  =  .9  and 

m  =  .001,  .1,  .2.  Varying  the  Xach  number  has  little  effect  save  for 
Y  in  the  region  of  -ll-h/HI. 


(i)  The  far-field  within  the  duct 


In  the  far-field  within  the  duct  a<x,u)  is  given  by  (A2.36)  of 
Appendix  A2  and  substitution  far  a  in  (2.29),  neglect  of  terms  of 
order  <oH/c)2ln(oH/c)  leads  to  the  following  expression  for  B(x, t); 


B(x,  t )  =  Ao(P(Y,  I»)-P<Y,  I- ) )  , 


(2.35) 


Ao  =  (1-V/Oc )  (Y/2x)  (x/2)*U*/B  , 


where 


P(t,x)  =  P» (t, x)-sgn(t)P- (t, x)  , 


P«(t,x>  =  q*  (t,  x)  +  JUdiq*  (t,x)Iln(H,  /8x(xa+t2>“)±g  (t,  x)3  )  , 


9- 


Figure  4 


B/4 m  of  (2. 33)  as  a  function  of  Odti/B.  h/H  1g  as 
Indicated  on  each  curve  and  (a)  JU  =  .  1;  (b)  .001. 
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Figure  5 


B/in  of  (2.33)  plotted  against  0«ltJ/B  for  h/H 
*«  =  .001,  .1,  .2  as  indicated  oi  each  curve. 
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q.(t.x)  =  I  ((xa+ta)“±x)/(xa+t3)l“  . 


(2.36a-c) 


(c.f.  Gradshteyn  and  Ryzhik  1980,  p484  and  p574).  g  Is  defined  by 
(2.34b).  Examination  of  (2.35)  and  (2.36)  shows  that  when  T  Is  large 
and  negative  (the  vortex  is  far  upstreaa)  the  field  in  the  duct  is 
small  with  B  a  Ao2/ I Yl *  while  when  the  vortex  is  far  downstreaa  of 

the  edge  B  *  -  Ao2(h/H)/ 1  Yl 3'2,  h<H  and  B  a  -  Ao2/IYI3'3,  h  >  H.  tfhen 

the  vortex  approaches  its  miniaum  distance  froa  the  plate, 

Ucttl/H  <<  t H-hi ,  we  find  for  K<  <<  1 

B/Ao  a  21  (1/1  1-h/HI  ••)  -  (1/  (l+h/H)**l  . 

and  B/Ao  has  a  aaxlaua.  This  can  be  seen  in  figure  6  which  shows 
B/Ao  as  a  function  of  IMtl/H  for  K,  -.1  and  various  values  of  h/B. 

As  for  the  radiated  far-field,  I  Bl  is  slightly  greater  when  the  vortex 
passes  above  rather  than  below  the  plate.  The  Influence  of  1U  on  B/Ao 
is  shown  in  figure  7  where  h/H=.9  and  1 U  =  . 01..1..2. 

In  figure  8  the  dependence  of  the  radiated  and  duct  far-field 
stagnation  enthalpies  on  Octtl/H  (respectively  B/Ae  of  (2.33)  and  B/Ao 
of  (2.35))  is  compared.  Two  values  of  h/H,  .5  and  .75,  are  shown  and 

*•=■!•  The  difference  in  the  two  is  aost  apparent  at  0«tt)/H  x  I 1-h/HI. 


82.6  Harmonic  Gust 

An  incompressible  velocity  perturbation  u.  =  (ui,u2)  of  a  convected 
two  dlaensional  gust  is 


IS 


Figure  6 


B/*o  of  <2.35)  as  a  function  of  O.tt  1  fB  with  X< 
values  of  h/H  as  indicated  on  each  curve. 
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Figure  8 


Caspar 1 sod  of  B /tm  of  <2. 33)  (solid  curves)  and  B/io  of 
(2.35)  (dashed  curves)  with  m  =  .l  and  h/H  -  .5,  .75  as 
indicated  on  each  curve. 
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I  =  <C<K>/K<>  (-Kz.Ki  >*zp(lK<  <Si-0«t>MK*h>,  (2.3?) 

•Aar*  C<K>  Is  a  function  of  the  wnve-nuaber  1  =  (X,,I2>.  Is  the  co¬ 
ordinate  systes  of  figure  1  the  vortlcity  of  the  gust  is  t|  (a.  Is  a 
uslt  vector  out  of  the  plane  of  the  paper)  asd 

f  =  (C<K>U.3*taa)/K,]exp(iK,  (i.-O.tHlbb)  ,  (2.38) 


'll* 


=  lC(i)<I,a+X2a>/t,l|  [exp(U,Xi+iIaX2)4  (x.-D.t-I.  )4  (x2-Xa>dX,dXa. 

(2. 39) 


The  far-fleld  stagnation  enthalpy  produced  by  the  Interaction  of  the 
gust  with  the  edge  (0,H>  of  the  plate  can  be  calculated  by  use  of  the 
results  of  S2.4  If  we  set 


h=Ia,  to~'*X,/0«t 


(2.40) 


and 


1  =  IC(I>(X,a*iaa>/X,lexp(lI,X,nXaIa>dI,dI2  .  (2.41) 


B(*>  Is  found  by  integration  of  the  stagnation  enthalpy  for  the  line 
vortex  case  of  82.4  over  the  variables  Xi,  Xa  as  in  (2.39).  Is  this 
way  we  obtain  with  (2.31) 


B(X)  =  C(l)(  1-T/0. 1  !P*(l-lsgn<Ki )  >a<x,  K,  Uc  )exp(-K,  0«t  tlbi/cJ  ) 
2*  11,1“ 


.  <ezp(lKa£)  -erp(-l  I,  IH>) . 


IX 


(2.42) 


■•suit  <2. 31),  used  to  derive  (2.45),  is  valid  for  short  vortex  waves. 


Thus  we  require  I l, HI  >>1  and  neglect  the  final  tern  la  the  curly 
brackets  of  (2.45). 

(!)  Outside  the  duct,  111 >>1/11.11 

from  (2.21)  with  the  asyaptotic  fora  of  the  Hankel  functions  and 
(2.42)  we  find 

B(j)  *  iC(I)U-¥/U.>  <  1 +cos8>  m,*sgn  <1.  )exp( IK2B-II1 0C  [  t)  )  ,  (2.43) 

2«(*l  xl  /■)*•(  ltiK.I<Hf(IiDc)/a] 

where  (tl  is  defined  by  (2.32)  with  to  =  0. 

(11)  Inside  the  duct,  l*im/II,II 

Fron  (2.45),  (42.36)  of  Appendix  42  we  find  for  the  far-field 

enthalpy 

B(*)  -  C(i)(l-¥/0c)<l-lsgn(4,))exp(i4aH-iKi0.lt)>  ,  (2.44) 

2«(IIi  I  H)*l  1  +  iK.LBf  (K.U<)/xl 

where  It)  is  defined  by  (2.32b)  with  to  *  0  and  6=0. 


42. V  Acoustic  energy  flux 


The  flux  of  acoustic  energy  (the  Blokhiatzev  definition,  c.f. 
Blokhiatsev  1946)  through  a  surface  S  is  Tl  say,  where 


32 


(2.45) 


TT  =  J\p*B  +  Op  •  B) .  a.  dS  . 


p,  II  =<U,0)  are  respectively  the  scan  density  and  velocity;  the 
perturbation  density,  velocity  end  stagnation  enthalpy  are  p* ,  i  and 
B.  a  is  a  unit  noraal  to  the  surface  S.  The  adiabatic  relation 
between  density  and  pressure  and  the  linearized  Bernoulli  equation  are 
used  in  (2.45)  to  give,  neglecting  IP  relative  to  unity. 


B(i  ♦  BB/c2) .  n  dS 


(1)  Line  vortex 


The  power  flux  through  the  surface  S*  of  figure  9,  a  circular  arc 
centred  on  the  origin  with  radius  S  =  )*)-•,  is  fl*.  say  where 


w 

n«  =  pj  B  iV„ 


+  R(cos8)B/c) I xl d8  , 


and  Bo  =  arsinlH/l si )-0.  v„  is  the  fluid  perturbation  velocity  in  the 
radial  direction  and  is  calculated  fros  Bernoulli's  equation; 


iv„/it  =  -\bA  I  xl 


(2.32)  and  expression  (2.33)  for  B  give, where  we  now  set  0.  =  D, 


iv„/it  =■  (l/c(l  +  *cas(8>))iB/4t  .  » 


Dslng  (2.49)  in  (2.47)  and  again  neglecting  IP  relative  to  unity,  we 
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Figure  9. 


TTr 


/ 


- - -Sr 

V- 


TT0 


'  r 


7777777  //'/,*  /  %  )}))}}'}}  FTTm 


Illustration  of  the  surfaces  in  the  far-field  through  which 
the  radiated  and  duct  power  fluxes  are  calculated. 
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TT*  »  (p/c)jBalxl  <18  i 

o 

and  with  (2.33)  this  becomes  with  0e  =  0, 

TT««/I  =  (3(l-»/0)alP/2)tF(0ttl/H1  1 1-h/HI  )-F(0ttI/H,  1+h/H))2  ,  (2.50) 

where 

I  =  p(y/2x)2/(H/0>  .  (2.51) 

The  8  dependence  of  (tl  has  been  neglected  in  obtaining  (2.50)  where 
now  (tl-  t-to-lxl/c.  The  flux  through  the  surface  S  at  ii  =  ®,  x2<H, 
(see  figure  9)  is  Flo  say,  and  fraa  (2.46),  (2.35)  we  obtain  for  IP<<1 

nc/I  =  («/2>(l-¥/U) 2H1  P(0(  tl  /H,  I  1-h/HI  )-P (Ut  tl /H,  1+h/H) 1 2,  (2.52) 

where  ttl  =  t-to-lxl/c  . 

In  figures  10,  11,  and  12  the  power  fluxes,  noraalised  by  I  of 
(2.51)  are  platted  in  dB  against  D[t)/H.  IDItl/HKl,  and  V/0  =  0.6. 

TTo  and  TT*  are  cowpared  in  figure  10(a-c)  where  *  =.l  and 
(a)  I 1-h/HI  =  .25;  (b)  .1;  (c)  .01.  Each  curve  has  a  singular  point 
because  TTw/I  =  0  where  B/in  of  (2.36)  is  zero  and  logio(TTw/I>+ 
Similarly  TTo/ I  =  0  when  B/Ao  of  (2.35)  is  zero.  Over  the  range  of 
01  tl  shown  in  the  figure  we  see  that  the  duct  power  flux  is  greater 
than  that  radiated  except  for  a  period  in  It]  >0  which  contains  the 
point  TTo  =  0  (e.g.  for  h/H  =1.25  in  figure  10(a),  this  range  is 
.31  <  OCtl/H  <  .75).  The  difference  in  the  results  for  the  vortex 
passing  above  and  below  the  plate  (if  there  is  any)  is  Indicated  in 
each  of  the  plots. 
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Figure  11.  lOlog.oOWf)  of  (2.50)  as  a  function  of  (MtJ/H  with 
V/0  =  .6,  K=.l  .  h/H  is  indicated  on  each  curve. 
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-  10  log)0  ( irR  /  N  ) 

- 10  io9l0rrr0/  n  j 


h/  H  =  0  »_  V/U  »  0  6 


Figure  12.  101ogio(T1n/I>  (solid  curves)  and  101ogio (FI  D/I>  (dashed 
curves)  plotted  against  Odtl/H  with  V/0  =  .6  and  ](  =  .01, 
.1,  .2,  as  Indicated  on  each  curve. 
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Figure  11  showB  TTn/I  for  different  values  of  h/H  and  again  M=.l. 

The  peak  of  each  curve  increases  as  the  miniaua  distance  of  the  vortex 
froa  the  plate  IH-hl  Is  decreased.  Though  the  functions  F  and  P 
(respectively  of  (2.34),  (2.36)),  appearing  in  the  radiated  and  duct 
power  fluxes  have  a  Mach  nusber  dependence  the  predoalnant  influence 
is  that  TTr/I  varies  as  IF*  whereas  TTo/I  is  proportional  to  X.  This 
is  illustrated  in  figure  12  where  TTr/I  (solid  lines)  andTTD/I 
(dashed  lines)  are  plotted  for  h/H  =  .9  and  X  =  .01,  .1,  . 2  as 
Indicated  on  each  curve. 

(11)  Haraonic  gust 

TT„,  TT0,  defined  in  the  preceding  section,  are  respectively  the 
power  radiated  to  the  far-field  outside  the  duct  and  through  the  duct. 
In  a  calculation  similar  to  that  described  between  (2.45)  and  (2.50) 
with  B  now  given  by  (2.43),  (2.44)  (additionally  taking  real  parts  of 
functions  and  averaging  over  a  wave  period  2*/K<U  )  the  following  are 
obtained: 


TTr  =  3pHI  C (K>  I  2(1-V/U)2X  ,  (2.53) 

16c*2 1  HiK,XHf(K,H)/*l2 


nD  =  pHIC(g) l2(l-V/B)2  ,  (2.54) 

2c*2K, HI 1+iK, XHf (X, 0) /*! 2 


and 


Tlr  =  TlR/fTo  =  3XIK,  IH/8  , 


(2.55) 


=  3lulB/8c  . 


where  u  =  1.0  Is  the  radian  gust  frequency. 


S3  Plane  wave  radiation  froa  the  duct 


S3.1  The  far-fleld  potential 

The  radiation  of  a  plane  wave,  f,exp(-lut),  froa  the  seal-infinite 
duct  foraed  by  the  plate  and  wall  is  now  discussed  (see  figure  13). 

The  incident  plane  wave  is  given  by 

ft  =  0oezp(-iukoX<  /  <1-10  )  ,  ko  =  o/c  ,  (3.1) 

and  propagates  upstreaa  froa  within  the  duct,  fa  is  a  constant  and 
flow  quantities,  such  as  X,  p  and  c  are  as  described  in  §2. 


The  total  perturbation  potential,  (fexp(-lut)  say,  is  the  solution  of 
the  convected  wave  equation  ((2.2)  with  the  right  hand  side  replaced 
by  zero)  which  satisfies  the  boundary  conditions 

r  V(o)exp(i)(X, )  ,  xi  >  0,  xa  =  H±0, 

=  < 

l  0  ,  Xa  =  0,  I X,  I  <  “  , 

and 

f  ^  constant. exp(ikolxl/(l+Xcos(9))/lxl“  ,  •  , 

0  =  arcos(xi  /  I  xl ) , 

for  x  outside  the  duct,  f  and  ^Axa  Bust  also  be  continuous  on 
X,  <  0,  Xa  =  B.  The  boundary  conditions  in  (3.2)  are  the  saae  as 
those  for  the  vortex  which  are  dlcussed  in  §2.2.  The  Prandtl-Glauert 


(3.2) 


(3.3) 


DISPLACEMENT  WAVE 


Figure  13.  todietlon  of  «  plane  «ve  fro.  the  duct 
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transform  of  (2.9)  is  employed  again  with  i  denoting  the  transform  of 

A J 

i  and  ft  that  of  ft  where 


Si  =  Soezpf-ikoXi )  . 


With  the  Green's  function  of  equation  (2.11)  we  find  on  application  of 
the  Divergence  theorem  to  the  fluid  in  the  region  V  sketched  in 
figure  3  that 


S(X>  = 


j<0< 


x;y>9s<y>  -  f(y)?G(x;y)  )■  a  ds  , 


where  S,  the  boundary  of  V  is  as  described  in  S2.3  and  sketched  in 
figure  3.  The  radiation  conditions  satisfied  by  G  and  S  ensure  that 
the  portion  of  the  surface  integral  over  Sr,  the  circular  part  of  S, 
is  zero  and  that  only  the  component  ft  of  0  contributes  to  the 
integral  over  the  surface  y,  =  «,  ly2l  <  H.  On  applying  the  condition 
(3.2)  we  obtain  with  (A2.4)  of  Appendix  A2 


t  (x)  =  -))(o)j  IG(x;y)2®*p(lxy> >  d7>  ~  0o2ikoKA(x. o>. 


The  notation  Ifl  is  explained  by  (2.6)  and  G(x;y>  -»  A(x,  u)exp(ikoyi ) , 
y i  -»  •»  defines  A,  which  is  calculated  in  Appendix  A2.  The  integral  on 
the  right  hand  side  of  (3.6)  arises  from  the  presence  of  displacement 
waves  on  the  plate  boundary  layers.  Ve  assume  that  since  the 
displacement  waves  must  decay  or  become  incoherent  as  they  propagate 
downstream  the  major  contribution  to  the  integral  comes  from  the 
region  y^  <<  1/y,  and  the  Integral  is  evaluated  by  use  of  the 


f  *T  ♦ 

approximation  (2.20a),  IG]_  =  atji,  which  Is  valid  for  I  yl  within  an 
acoustic  wavelength  of  the  plate  edge,  Ixl  in  the  far-field.  On  the 
assumption  of  short  displacement  waves,  x&  >>  1,  a  further 
approximation  for  [£)  given  by  (41.11)  can  be  used  in  (3.6)  and  with 
result  (2.24)  we  find  on  reverting  to  the  actual  potential  f 
(c.f.  (2.9)) 


A(X>  =  -exp(-ikoKxi ) 


a(x,u)|/(<j)  (H/l  xl  )*(l-i)  +  0o2ikoHA (x, u) 
l  IXI 


(3.7) 


J)(o)  is  determined  by  the  Kutta  condition  in  Appendix  A4,  (44.16)  and 
a,  A  are  given  respectively  in  equations  (A2. 12) , (A2 . 13)  of 
Appendix  A 2  far  x  in  the  far-field  outside  the  duct  and  by  (A2.17), 
(A2.14)  when  x  is  far  downstream  within  the  duct. 


Substitution  for  V(o)  gives 


*<X>  =  -koH^oexp(-ikoMxi )  (  (( l-R+ll{l+B)la(x,o)/x)  +  2iACx,o))  . 


(3 .  B) 


6  is  a  reflection  coefficient  when  0  =  the  perturbation  potential 
in  the  absence  of  displacement  waves  and  is  defined  by 


fa(x)  -*  #1<X>  +  #r(x)  ,  X.  -*  •,  I  Xa  I  <  H, 


(3.9) 


t r  (x>  =  AoBexp(ikoXi/(l+K))  . 


B  is  determined  from  (3.7)  by  setting  |)(u)  =  0  and  use  of  results 
1*2.14),  (A2.15)  of  Appendix  A2.  Then  comparison  with  (3.9)  gives 


1  -iuHf(u)/xc 


(3. 10) 


1  +luHf(o)/xc 


where  f  (<j>  is  defined  la  (2.21b). 


(1)  The  radiated  field 


As  in  (2.7)  we  set  x  =  o/V.  Results  <A2.8a,b>,  (A2.12),  (12.13)  of 
Appendix  A2,  in  which  the  asyaptotic  fora  of  the  Hankel  functions 
(Abraaowltz  and  Stegun  1964,  p364)  are  substituted,  give  with  (3.8) 
and  (3.10)  for  olxl/c  -*  • 

f  (x)exp(-iut) 


->  p<,(l-l>  (oH/xc)“ 


(I  xl  /H)' 


1- (ioHf (w)cos0/xc) (l+cos8 ) 


( 1+luHf (u)/xcl 


exp(-lu(  t) ) , 


(3. 11) 


c 


<V/U)H 


fl+iuHKf (u)/xc 


1+iuHf (u)/xc 


(3. 12) 


where  it)  =  t- I xl / (c ( l+McosO) ) ,  correct  to  neglect  of  teras  of  order 
JP  relative  to  unity.  The  far-fleld  pressure,  pexp(-iut),  is  found 
froa  the  linearized  Bernoulli  equation  which  gives,  again  neglecting 
IP  relative  to  unity,  p/p  *  lu(l-Rcosfi)p. 

The  intensity  of  the  radiated  sound,  averaged  over  a  wave  period 
2x/o,  is  lpl*/2pc  and  froa  (3.11)  we  find 


pi*  = 


Ip.  I 
2pc 


2a-I)*koH  ft.  ,  olxl/c  >>  1  , 


<3. 13) 


2pc 


2pc  J 


<xl*l/B) 


where  Du  is  a  directivity  factor  defined  by 


ft.  =  <l-Ilcos8>*|  1-  (ioHf  <o)  <cos8) /xc)+c  <l+cos8) 


1  +  ioBf  <u) /xc 


<3. 1«) 


and 


I p, I */2pc  =  pftil pol / <l-M))a/2c 


(3. 15) 


is  the  intensity  of  the  incident  wave  f  ,  of  (3.1).  In  the  long 
wavelength  limit  we  find,  neglecting  of  terms  of  order  (oB/c)* 
relative  to  unity 


ft.  =  < l-Ecos8)*  <  1- ( 1 +coe8  >[  k«B- ( V/U ) B(2-hoH (3+coe8) ) 1 )  ,  (3.16) 


oH/c  <<1. 


The  term  involving  V/0  in  (3.15)  is  the  change  in  the  far-field 
intensity  from  that  of  the  'no  displacement  wave  case'  and  indicates 
that  for  long  sound  wavelength  the  intensity  is  increased  by  an  amount 
which  is  proportional  to  VH/U. 

The  directivity  factor  ft.  is  plotted  against  8*  in  figure  14  for 
various  values  of  koH,  V/U  =  .5  and  (a)  ■  =  .1;  (b)  ■-  01.  The 
maximum  is  at  8  =  ISO*,  upstream  of  the  duct  on  the  well  at  ia  =  0, 
when  ft.  =  (1*I)2.  The  minium  attained  by  ft.  depends  on  koH, 
decreases  as  koB  increases  and  occurs  for  koB  I  1  at  0*.  The  results 
presented  here  are  strictly  valid  only  for  low  frequencies,  certainly 
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FIG  14  (a) 


Figure  14.  The  directivity  factor  of  tha  far-flald  pressure,  U.  of 
(3. 18),  as  a  function  of  8,  with  f/0  =  .8  and  (a>  M=.l; 
<h)  .01.  The  value  of  koH  la  indicated  on  each  curve. 
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below  the  cut-off  of  transverse  duct  aodes,  l.e  koH  <  x/2,  but  results 
are  shown  for  greater  values  of  koH.  The  figure  reveals  that  Dr  can 
become  very  small  for  values  of  koH  >  1  and  6  close  to  0*.  Figure 
14(b)  shows  that  between  the  upstream  and  downstream  directions  there 
is  a  difference  of  4  dB  at  koH  =  .5. 


<ii )  The  duct  field 


The  far-field  potential  within  the  duct  can  be  expressed  as 

rf(X>  =  4t(x>  +  4=.(x)  +  4*>(x>  .  (3.17) 

where  4t>(x)  is  a  reflected  wave  arising  from  the  presence  of 
displacement  waves  on  the  plate  and  4a(x>  is  the  reflected  field  in 
their  absence  given  by  (3.9),  (3.10).  The  total  reflected  field, 

<X)+4*»<X)  can  be  written  as  a  single  reflected  wave  with  reflection 
coefficient  R; 

(x>+4s<x>  =  4oR  exp(ikoxl/U+K))  ,  (3.18) 

and  from  (3. 8)- (3. 10) 

B  =  -  <  1  -  (ihoflf  (u)/x)  +  2e  >.  (3.19) 

(  1  +  ikoHf  (o)/x  > 


In  the  absence  of  displacement  waves  e  =  0  and  then  in  the  long 
wavelength  limit  R  ^  -1.  In  terms  of  pressure  amplitudes  with  p,  that 
of  the  incident  wave  and  pr  that  of  the  reflected  field  of  (3.18)  we 

find 


S’* 


T 


ip.i/ip.i  =  <  <i-m  /  >  i  ri 


(3.20) 


Experiments  have  been  conducted  at  the  Institute  of  Sound  and 
Vibration  Besearch  of  Southampton  university  to  measure  the  field 
reflected  at  the  mouth  of  an  open  pipe  which  has  a  low  mean  flow  Kach 
number  Inflow  and  white  noise  excitation  within  the  pipe,  downstream 
of  the  mouth  (Davies  1967).  Comparison  with  the  present  theory  is 
possible  only  on  a  qualitative  level  since  in  the  experiments  the 
geometry  is  that  of  a  circular  pipe,  rather  than  a  two-dimensional 
duct;  the  inflow  is  produced  by  suction  within  the  pipe  and  is 
consequently  non-uniform  around  the  mouth.  Also  the  theory  assumes 
that  the  flow  does  not  separate  at  the  duct  mouth  whereas  for  a  sharp 
pipe  entry  the  experimental  flow  consisted  of  a  shear  layer  of 
thickness  around  ,2a  (a  is  the  pipe  radius).  A  bell  mouthed  pipe 
entry,  also  examined  in  the  experiments,  had  smooth  entry  flow  but  its 
pipe  wall  (of  radius  of  order  a  at  the  mouth)  does  not  resemble  the 
thin  wall  discussed  here.  levertheless  in  the  latter  case  the 
measurements  of  the  ratio  of  equation  (3.20),  normalized  by  its  no¬ 
flow  value  are  curve  fitted  by  ( (1-K) / (1+K) 3 ■ 9 •  which  Davies  points 
out  is  £  eater  than  the  value  I  (1-JC)/ (1+K)  1  which  a  potential  flow 
solution  would  predict.  The  presence  of  coupled  instability  waves  in 
the  pipe  boundary  layers  is  suggested  as  the  cause  of  the  enhanced 
reflected  wave  amplitude.  This  explanation  is  in  accord  with  the 
present  analysis  which  from  (3.12),  (3.19),  (3.20)  gives  the  ratio  as 

( 1— H) 1 1+2 (V '0)  10  =  (  U-10/(1+X)J<’-'",J’  , 


(1+*) 
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(R-  and  (koH)2  are  neglected  relative  to  unity).  This  prediction  of 
the  ratio  is,  for  7/0  =  .6,  such  greater  than  that  of  the  experiaent 
but  shows  that  flow/ acoustic  interaction  at  the  duct  south  can  explain 
the  Increase  in  the  reflected  field. 

3.2  Acoustic  energy 

Equation  (2.45)  gives  the  instantaneous  acoustic  energy  flux  through 
a  surface  S.  Here  we  discuss  the  tine-averaged,  i.e.  averaged  over  a 
wave  period  2x/o,  power  flux  away  froa  the  plate  leading  edge.  The 
flow  of  acoustic  energy  through  the  surface  S,  sketched  in  figure  9 
and  defined  by  x>  =  L  >>  c/o,  Ixzl  <  H,  is  TTn.  The  perturbation 
field  at  S  is  given  by  equations  <3. 17)- (3. 19) .  Restoring  the  tiae 
factor  exp(-lot),  taking  real  parts  of  functions  and  averaging,  yields 

TTo  =  -  T1  iU-IRI*>.  (3.21) 

where  TTi  is  the  power  flux  of  the  incident  wave  y,  of  (3.1)  and 

TT i  =  pol«ol2koH/2  .  (3.22) 

Calculation  of  IRI2  froa  (3.19)  correct  to  second  order  in  koH  shows 
that 

Tlo/fli  =  -2koH(  1-koH  -  2(V/0)M  l-2koHl/koH  )  ,  koH  <<  1.  (3.23) 

The  tiae-averaged  flux  through  Sr  of  figure  9  (a  circular  arc  centred 
on  the  origin  with  kolxl  )>  1)  TTr  say,  is  found  in  a  slailar  Banner 
using  (3.11)  and  (3.12); 


S3 


Tu/TTi  =  2koH  (  ll+£la  +l£-ikoHf (u)/xl2/2  )  , 


(3.24) 


I  1  +  ikoHf (o)/xl 2 

where  £  is  defined  by  (3.12)  and  f  by  (2.21b).  In  the  limit  koH  -»  0, 
neglecting  IP  relative  to  unity,  this  becomes 

rr„/n  4  -»  2koH  <  1-koH  +2X(¥/U )  )  .  (3.25) 


In  the  long  wavelength  Unit  the  net  flux  of  acoustic  energy  away  from 
the  plate  edge  is  fron  (3.23)  and  (3.25) 

TT*  +  Tlo  =  TT  ,4H(V/U)  (  1-koH  >  ,  koH  <<  1.  (3.26) 

Acoustic  energy  is  not  conserved  due  to  the  presence  of  displaceaent 
waves  on  the  plate  boundary  layers  and  (3. 26)  shows  that  for  low 
frequencies  the  total  flux  away  fron  the  plate  is  positive;  sound 
energy  is  produced  by  the  flow  acoustic/ interaction  at  the  plate  edge. 


In  figure  15  the  radiated  and  duct  energy  fluxes  normalized  by  Tl  i 
and  respectively  of  (3.21),  (3.24),  are  plotted  in  dB  against  koH  for 
V/U  =  .6.  Also  shown  is  a  transmission  factor  T1  t  =  FTr/HTdI,  which 
is  the  radiated  power  as  a  fraction  of  the  duct  power  flux.  Though 
the  results  are  valid  only  for  koH  <<  x/2  a  much  greater  range  is 
shown  in  figure  15(a)  where  X  =  .  1.  A  more  detailed  picture  of  these 
results  for  koH  <  1  is  given  in  figure  15(b).  Singular  points  in  the 
curves  for  TTo  and  TTt  are  evident  at  koH  -  .1.  This  is  where 
I8I2  -  1  is  (3.21)  and  TTo  =  0.  For  koH  <  .1  the  net  duct  energy 
flux  is  positive  in  the  direction  of  the  reflected  field,  i.e.  in  the 
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FIG  15  (  a ) 

Figure  15.  The  normalized  radiated  and  duct  power  fluxes,  Tin 
Flo  of  (3.21),  (3.24)  respective! j,  plotted  in  dB 
against  koH  with  V/D  =  .6  and  (a>  X  =  .  1,  koH  i  4. 
(b)  X  =  .1,  koH  (  1.;  (c)  X  =  .01,  koH  (  1. 


-SS- 


and 


M  (  UH  A  >47411 


direction  of  the  inflow,  while  for  koH  >  .1  the  net  flu*  is  in  the 
incident  wave  direction.  The  corresponding  results  for  E  =  .01  are 
shown  in  figure  15(c).  Figure  16  shows  that  the  net  flu*  of  acoustic 
energy  away  fro*  the  plate  edge,  T7„  +  TId,  normalized  by  17,  is 
positive  as  Indicated  by  (3.26). 


FIG-  16 

Figure  16.  The  net  normalized  power  flu*  away  from  the  duct  mouth 
(rU+rTo)/n,.  as  a  function  of  koH  with  V/U  =  .6  and 
E=  .1,  . 01  as  indicated  on  each  curve. 


A 


Contusions 


The  production  of  sound  when  a  weak  line  vortex  Is  convected  into  a 
seal-infinite  two-dimensional  duct  by  a  low  Mach  nunber  nean  flow  has 
been  examined  by  use  of  a  low-frequency  Green's  function.  Howe's 
theory  of  displaceaent  thickness  fluctuations  was  applied  to  model  the 
effect  on  the  radiated  sound  field  of  disturbances  which  are  produced 
by  the  interaction  of  the  vortex  with  the  leading  edge  of  the  duct 
wall  and  which  propagate  in  the  thin  boundary  layers  emanating  fron 
the  edge.  The  strength  of  the  disturbances  is  fixed  by  the  Kutta 
condition  applied  at  the  leading  edge  and  It  is  found  that  the 
presence  of  displaceaent  waves  reduces  the  radiated  pressure  field  by 
a  factor  of  about  4/10  compared  to  the  case  in  which  their  presence  is 
ignored. 


The  duct  width,  H,  is  assuaed  small  relative  to  the  characteristic 
acoustic  wavelength  d/H,  where  d  is  a  varticlty  length  scale  and  H  is 
the  flow  Mach  number.  Analytical  expressions  are  obtained  for  the 
far-field,  both  within  and  outside  the  duct.  From  this  solution  we 
find  the  radiated  field  produced  when  a  frozen  harmonic  gust  is 
convected  by  the  duct  Inflow. 

The  instantaneous  power  fluxes  of  the  duct  and  radiated  fields, 
produced  by  the  interaction  of  the  line  vortex  with  the  duct  mouth, 
are  compared.  Ve  find  that  the  former  is  of  order  Y*M  whereas  the 
radiated  power  flux,  an  order  of  magnitude  smaller,  is  proportional  to 
0f/2x  is  the  vortex  strength). 

The  radiation  of  a  plane  wave  from  within  the  duct  has  also  been 
examined.  The  net  flux  of  acoustic  energy  away  from  the  duct  mouth  is 
found  to  be  positive  due  to  a  production  of  acoustic  energy  by  the 
displacement  waves. 
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Figure  captions 


Figure  1. 
Figure  2. 
Figure  3. 


Figure  4. 


Figure  5. 


Figure  6. 


Figure  7. 


Figure  8. 


Figure  9. 


Figure  10. 


Schematic  illustration  of  convection  of  a  line  vortex  by  a 
unifora  aean  flow  into  a  seal-infinite  duct. 

Sketch  of  the  idealized  boundary  layer  aodel  with  a 
displaceaent  wave. 

Sketch  of  the  control  surfaces  of  equation  (2.14)  and  the 
vortex  wave. 

B/An  of  (2.33)  as  a  function  of  Deltl/H.  h/H  is  as 
indicated  on  each  curve  and  (a)  JU  =  .1;  (b)  .001. 

B/A*>  of  (2.33)  plotted  against  Ueltl/H  for  h/H  =.9  and 
Me  =  .001,  .1,  .2  as  indicated  on  each  curve. 

B/Ao  of  (2.35)  as  a  function  of  U.ttl/H  with  Me  =  .  1  and 
values  of  h/H  as  indicated  on  each  curve. 

B/Ao  of  (2. 35)  plotted  against  IMtl/H  for  h/H  =  .9  and 

Me  =  .01,  .1  ,  .2. 

Conparison  of  B/A«  of  (2.33)  (solid  curves)  and  B/Ao  of 
(2.35)  (dashed  curves)  with  Me=.l  and  h/H  =  .5,  .75  as 
indicated  on  each  curve. 

Illustration  of  the  surfaces  in  the  far-field  through  which 
the  radiated  and  duct  power  fluxes  are  calculated. 

The  normalised  power  fluxes  Tit./ 1,  TTo/I  of  (2.50)  and 
(2.52)  respectively,  plotted  in  dB  against  Odtl/H  with 
V/U  =  .6  and  M  =  .1:  (a)  I  1-h/HI  =.  25;  (b).l; 


U 


Figure  11. 


101og,om*/I>  of  (2.50)  as  a  function  of  O.ttJ/H  with 


V/0  =  .6,  M= . 1  .  h/H  is  Indicated  on  each  curve. 

Figure  12.  101og,o(nR/I)  (solid  curves)  and  101ogto(T1o/l)  (dashed 

curves)  plotted  against  0c(t]/B  with  V/D  =  .6  and  X  =  .01, 
.1,  .2,  as  indicated  on  each  curve. 

Figure  13.  Radiation  of  a  plane  wave  from  the  duct. 

Figure  14.  The  directivity  factor  of  the  far-field  pressure,  Dr  of 
(3.16),  as  a  function  of  6,  with  V/0  =  .6  and  (a)  H=.l: 

(b)  .01.  The  value  of  koH  is  indicated  on  each  curve. 


Figure  15.  The  normalized  radiated  and  duct  power  fluxes,  T1r  and 
TTo  of  (3.21),  (3.24)  respectively,  plotted  in  dB 
against  with  V/0  =  .6  and  (a)  X  =  .1,  koH  (  4.; 

(b)  X  =  .1,  koH  <  1.;  (c)  X  =  .01,  koH  <  1. 


Figure  16.  The  net  nornallzed  power  flux  away  from  the  duct  mouth 
CTTR+rTol/TK,  as  a  function  of  koH  with  V/0  =  .6  and 
X=  .1,  .01  as  indicated  on  each  curve. 


Figure  17.  The  transformation  Z  of  equation  (Al.l). 


Figure  18.  Sketch  of  the  regions  1-3  around  the  duct. 


Figure  19.  Sketch  of  the  flow  of  equation  (A3. 4) 


The  Green's  function  of  (2.11)  of  the  main  text,  G(x;y)  Is 
calculated  In  Appendix  A2  for  x  in  the  far-field.  Then  when  y  is  in 
the  near-field  of  the  duct  mouth  (BB'  of  figure  17)  G<x;y> 
approxiaates  the  potential  of  an  incompressible  flow  into  the  duct. 
Thus  the  potential  and  stream  function  of  the  potential  flow  into 
the  duct  which  has  a  uniform  velocity  of  unity  in  the  positive  x< 
direction  at  Xi  =  •,  lx2l  <  H,  are  required  in  the  calculation  of  the 
far-field  stagnation  enthalpy  and  are  discussed  below. 

The  mapping  of  the  z  plane  (z^xi+jx^  and  j  is  the  complex  imaginary 
number  (-1)“)  consisting  of  two  semi-infinite  planes  defined  by  x2  - 
±H,  x.  >  0,  in  the  (xt.Xa)  co-ordinate  system  in  figure  17  to  the 
upper  half  of  the  Z  plane,  also  sketched  in  the  figure,  is  found  from 
the  Schwarz-Christof fel  transformation  (Kilne  Thompson  1968, §10. 2); 

dz/dZ  =  KiZ^-ll/Z  ,  (A1 . 1 ) 

which  gives  on  integration  <with  the  complex  constant  K  and  that  of 
the  integration  chosen  so  that  the  points  <0,*H)  are  mapped  to  (±1,0) 
in  the  Z  plane) 

z  -  (H/xM  Z 2  -1  -ln(Z2>  +jx  )  .  (A1.2) 

A  uniform  flow  of  velocity  1  at  ii-*»,  lx*l<  H,  is  equivalent  to  a  sink 
at  the  origin  of  input  2H  taken  over  an  angle  x  and  therefore  of 
strength  2H/x.  The  complex  potential  of  this  flow  is  w  and 


Figure  17.  The  transformation  Z  of  equation  (Al.l). 


Figure  IB.  Sketch  of  the  regions  1-3  around  the  duct. 
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w  =  -12H/x)ln(Z>  +JH 


(A1.3) 


and  re-arrangement  gives 

Z  =  +jexp<-xw/2H) . 

Substitution  for  Z  in  (A1.2)  yields 

z  =  {H/*X-1  +  exp<-xw/H>  +  wit/H)  . 

The  fluid  velocity  u  =  (u,v>  is  found  from  u-jv  =  dw/dz  and 

dw/dz  =  1/(1  +exp(-*w/H>)  . 

Equating  real  and  imaginary  parts  in  <A1.5)  with  w  =  we 

x,  =  <H/x)  f-1  +  -  cos(*J/H)exp(-xJ/H)  }  , 

x2  =  <H/*X  (J'v/H)  +  sin<xJ/H)exp<-x£/H>  )  . 

Appr  axi  i  ions_J  or 

On  I,  !  0,  i;  =  ±H,  is  constant  and  ^  =  ±H.  As  Xi->°>  with 
lx2l  <  H  examination  of  (A1.7)  shows  that  £  ->  Xi  .  As  1*1  -* 
the  duct  region  of  figure  A1  we  find 

J  h  -  <H/x)  ln(l*lx/B)  . 

Expansion  of  the  exponential  factors  in  (A1.7)  for  e  <<  H/n 
x2  =  H  =  'J  leads  to 


bb 


(A1 . 4 ) 


(A1.5) 


(A1.6) 


find 


(A1.7) 


®  outside 


<A1. 8) 


with 


Zi 


=  (H/2x>  (j^/H>2{  1+  0(xJ/H>> 


a  <<H/x  . 


(A1.9) 


Thus 

£  =  -(±(2Hx,/x)“)  ,  j;  =  H±0,  0  <  x,  <<  H/x  , 

where  the  upper/lower  signs  correspond.  The  disconinuity  in  0 
the  surface  xs  =  H,  xi  >  0,  is  approximated  by 

[|t.  =  2  (2Uxi  /»)**  ,  x,  <<  H/x  . 


Appendix  A2.  Green's  function  calculation 


Define  G(j(;y)  as  the  solution  of 


-  i2  + 

.  iy,2 


u2 

c2 


G(x;y) 


which  satisfies  the  boundary  conditions 


dG(x;y)/^y2 


0.  yz  =  0,  ly,l  <  ®, 

0,  ya  =  H,  yi  >  0  , 


and  the  radiation  condition  that  for  toy/cl  >>  0(1),  G(j;y)exp( 
should  have  inconing  wave  properties  with 


(Al. 10) 


across 


(Al. 11) 


(A 2. 1) 


(A2.2) 


(A2.3) 


-iut) 


G(*;y)  -i  Aexp(ioy,/c)  ,  y,  ^  «,  ly2l  <  H. 


(A2. 4) 


In  y  space  C(x;y)  represents  a  source  with  position  x  and  will  be 
calculated  for  Ijl  in  the  far-field.  The  boundary  condition  (A2. 2) 
will  be  satisfied  by  placing  an  iaage  of  the  seni-inf inlte  plane  at  y2 
=  -H,  yi  >  0,  and  an  iaage  source  at  x*  =  (11,-12)  and  requiring 
instead  that  ^G<z;  y) /^y2  =  0  on  the  iaage  plane.  The  flow-field 
around  the  duct  foraed  by  the  two  planes  is  divided  into  three  regions 
<see  figure  18).  In  region  1,  olyl/c  <<  0(1),  while  uy,/c  >  0(1), 
ly*l  <  H  in  region  2  and  in  region  3  luy/cl  >  0(1).  We  consider  both 
the  case  of  x  in  the  far-field  outside  the  duct  and  within  the  duct. 

(i)  x  in  the  far-field  outside  the  duct 

Region  3  :  The  function  G(x;y)  will  take  the  fora 

C(x;y)  -1  G,(*;y)  +  G,(x;y) 

where  for  I  ox/cl  >>1 

(■)  (u 

G.(x;y>  =  (1/1)1  Ko(ulx  j[l/c)  +  Ho<<jI  x*-yl  /c> 

and  (A2.6) 

(.') 

Gs(x;y)  =  P  Ho(uiyi/c)  . 


U 

Ho(x)  is  a  Hankel  function  of  the  first  kind  and  p  is  an  unknown 
function  of  x  The  function  Gt(x;  y)  is  the  solution  of  (A2.1)  when 
only  a  plane  at  y*  =  0  is  present  while  Gy<x;y)  arises  froa  the 
presence  of  the  plates  and  is  the  far-field  produced  by  the 
scattering  of  Gt<x;y)  by  the  plate  edges.  As  luy/cl  -»  0, 


G,  (x;  y>  -»  a  +  byi  , 

and  (42.7) 

(^(Xjy)  -*  <2p/*){  iln(olyl/c>  +  * < . 5772-iln2) /2)  , 

where 

a  =  <i/2)Uo(ul  gl  /c) ,  b  =  (ioxi /2cl xl >H, (ol xl /c)  ,  (42.8a,b> 

(c. f .  Abraaowitz  and  Stegun  1964,  p360>. 

Region  1:  Here  loy/cl<<l  and  the  wave  equation  (42. 1)  is  approximated 
by  the  Helmhlotz  equation; 

{V  +  V  'l  G(x;y)  =  0,  I  xl  >>  lyl. 

I  ^y.* 

G(X;y>  will  be  of  the  form 

G(x;y)  =  a  +  byi  +  aa(y)  +  r,  (42.9) 

where  a,  tr  are  as  yet  undetermined  functions  of  x  and  J(y>  is  the 
potential  of  an  Incompressible  flow  which  satisfies  ^Ay2  =  0  on 
y^-±H,  yi  >0  and  as  yi  -»  •»,  ly2l  <  H,  a  -*  y, .  £  is  calculated  in 

4ppendix  41  (equation  (41.7).  Hatchng  the  far-field  fora  of  (42.9) 
with  the  'outer  field'  of  (42.6)  gives  with  (42.7)  and  result  (41.8) 
of  4ppendix  41 

2p(iln(ol yl /c)+x/2  +i . 5772-lln2)  -  a(-Hln(ol yt /c) +Hln(oH/xc) ) +  xv  . 


(42.  10) 


Region  2:  G(*;jr)  takes  the  form  given  in  (A2.4)  and  another  equation 
Is  obtained  by  Batching  the  near-field  form  of  that  equation  to 
G<x;y>  of  region  1  given  by  02.9)  in  the  llait  uyi/c  ->  «.  In  this 

May  we  obtain 


a  +  byi  +  oyi  +  r  =  A(l+iuy,/c)  .  (A2.ll) 

The  equations  (A2.10),  (A2.ll)  determine  a,  p,  r,  A  in  teras  of  the 
functions  a  and  b  given  in  02.6).  Ve  find 

a  =  (-bi  iua/c) /  (It  1  (oH/itc)  f  (u) )  ,  (A2.  12) 

and 

P  =  ioH/2  ,  r  -  -aUf  (u)  /x  ,  A  =  (a  +  b)/iiu/c>.  (A2. 13) 


where 

f <«j>  =  . 5V?2  +  ln(<jH/2*c>  -  ix/2. 

<ii)  x  in  the  far-field  within  the  duct 

In  a  calculation  siallar  to  that  outlined  in  <i)  above  we  obtain  for  A 
of  02. t) 


02 .  14) 


CA2. 15a, b) 


In  region  1  we  find  that 

G<x;y>  *  ♦  r<x>, 


A  =  Q  +  lexp<-luxt /c)/ (2uH/c) 

and 

A 

Q  -  -Rexpde'Xi /c)  ,  R  =  -  {1-ioHf  (u) /xc)  . 


2loH/c  (1  +  iuHf(u)/xc) 


?0 


<A2. 16) 


where  now 


a  =  <iuB/c)  ♦  eipduii  /c)/2H  =  1.  exp ( luxi  /c>  ,  (12.17) 

S<1(  loHf  (u)/xc) 

and  r(x>  is  a  function  of  x  which  is  not  needed  in  the  calculation  of 
the  wain  text. 


Appendix  13- 


the  boundary  layer  displaceaent 


waves 


Define  G(x;y>>  as  the  solution  of 


t  i2  +  V2  - 
I  ^x> 2  ixa2 


r  i  u 

l  c 

+  KA  1 

*  ' 

6  x ,  j 

0  , 


<13. 1) 


which  satisfies 


f  i(ii'fi),  Xa  =  H±0,  ii  >  0, 

ic/i»2  =  < 

.0,  xa  =  0,  I  It  I  <  ■  . 


(13.2) 


For  ulxl/c  <<  1  and  neglecting  convection  by  the  aean  flow  since  M<<1 
(13.1)  becoaes 


A2  t  ^2 
^Xi2  ixa2 


G<x;y> >  =  o  ,  uixi/c  <<  l. 


(13.3) 


1\ 


If  also  oly>l/c  <<  1,  (5(x;y)  represents  a  line  source  of  strength 
±l/%  respectively  on  the  upper/lower  side  of  a  rigid  plate  with  a 
rigid  wall  a  distance  B  below  (see  figure  19)  in  an  incompressible 
fluid.  In  the  Z  plane  (c.f.  Appendix  Al,  (A1.2>>  we  find  that 

G(x;y)  =  Ke(Va(x;y>)  +  a  ,  ulxl/c,  ulyil/c  <<  1  ,  (A3. 4) 

where 

Vo  =  llnfZ2-Z»2)  +  Xln(Z2)  .  (A3. 5) 

«  Z2-Z_2j  x 

1  and  a  are  constants.  The  term  in  X  accounts  for  flow  at  Xi  =  *, 

IXil  <  H  and  Z*  are  the  images  of  the  points  yi+jH±0  in  the  Z  plane. 


At  the  plate  leading  edge,  x  =  (0,H),  Z  =  1.  tfe  find  from  (A3. 5) 


dVa 

=  2 

i 

l 

,  ol yi 1 /c  <<1  . 

dZ  Z=1 

X 

(1-Z»2) 

(1-Z-2) 

(A3. 6) 

For  values  within  a  few  duct  widths  of  the  duct  mouth  x  is  small 
and  has  been  neglected  in  (A3. 6)  where  the  remaining  terms  are  (c.f 
Appendix  Al)  of  order  (xyi/H)-**  for  yi  <<  H.  The  component  of  the 
perturbation  potential  due  to  the  presence  of  displacement  waves  on 
the  boundary  layers  of  the  plate  at  Xa  =  tt±0,  Xi  >  0  can,  in  view  of 
(A3.1),  (A3. 2)  and  equation  (2.8)  of  the  main  text  be  expressed  as  a 
complex  potential,  Vb  say,  in  the  near-field  of  the  plate  leading  edge 


with 


The  perturbation  potential  due  to  the  vorticity  wave  u>  (<2.4b>  of  the 
aain  text)  will,  in  the  vicinity  of  the  plate  edge  <0,H),  approx  1  mate 
that  of  an  incompressible  flow  which  has  a  complex  potential  Vd  where 


r 

Vd  --  -J  J  exp<io(to-t) )  I 

fm 

[Z-Zo 

2x(2x0.)  J 

tz-Zo* 

explloy  i  /U.  )dyi 


-J  t  exp<io<t0“t) )  j  In 
2*  (2x0. ) 


Z+Zo 


ZtZo 


exp(iuyi /U. )dy i  .  (44. 1) 


Zc.  is  the  image  of  the  point  yi+Jb  under  the  mapping  (41.2)  of 
4ppendix  41.  The  second  integral  on  the  right  hand  side  of  (44.1) 
accounts  for  the  image  of  the  vortex  wave  in  the  wall  at  i2  0.  The 
total  potential  in  the  vicinity  of  the  edge  <0,H)  has  complex 
potential  V  say,  where  from  Appendix  43,  (43.9)  and  (43.4) 


V  =  V„  +  Vs.  <44. 2) 

Vs  arises  from  the  boundary  layer  displacement  waves  while  Vd  is  due 
directly  to  the  vorticity  wave.  The  fluid  velocity  at  the  leading 
edge  of  the  plate  is 


dV 

d  (Va  ♦  Vs) 

dZ 

— 

— 

— 

dz 

JM  < 

dZ 

<±Z 

-  I  <x/2H> (d <Vd  1  Vs) /dZ> / (Z2-l) 1  I Z= 1  ,  (44.3) 


< 


Cc.f.  (A1.2)  of  Appendix  11).  For  the  fluid  velocity  to  remain  iinite 
at  the  edge  we  require  that 

d(  +  Wt,  )/d2  =  0  ,  Z  =  1.  <A4.4> 


Froa  equations  (A3. 6),  (A3. 9)  of  Appendix  A3  we  obtain  for  Imty)  >  0 


(dVe/dZ) I Z  =  1  =  2(^(o)/*)exp<iot) 


exp ( i xy i >d£ 


<A4.b> 


where  the  integration  variable  is  the  potential  0  and  (A1.4),  (A1.6) 

01  Appendix  A1  have  been  used.  With  (A1.7)  of  Appendix  A1  and 
Gradshteyn  and  Ryzhik  (1980,  p308)  we  find 


<dV„/dZ) I Z-l  -  -2  (R  (olH/xd  expi  L  <- lyH/x) i  ,  <A4.6> 


where 


L(x)  -  lnT(x)  -xlnx  t  x  . 


(A4. 7) 


T(x>  is  the  gamma  function  (Gradshteyn  and  Ryzhik  1980,  p933). 


Differentiation  of  (A4.1)  gives 


dWd 

dZ 


=  <2/xD« )  {y/2x>exp<lo(to  t) , 
Z  1 


<A4. 8) 


where 


K- 


expdoy  1  /0.  )  Iaf  1/ (l-Zo*>  1  dy. 


<44. 9) 


3S 


lm(f)  denotes  the  imaginary  part  of  f  and  Zo  Is  the  image  of  yi«Jh  in 
the  Z  plane.  Condition  (14.4)  becomes  with  (14.?),  (14.8) 


^(o) 


(t/2HU,  )expl  iot0-L(-ixH/«)l, 


(14. 10) 


lELjevaluAlt 


& 


■ate  that  (c.f.  Appendix  Al,  (A1.4),  (A1.6>) 


Jl 


T(o)  +  T*(-o)  , 


(14. 11) 


where 


T  (u) 


-iexp(oh/U, ) 

2  ,J 


dw  e*p( iuz/Ut )dz 
dz 


(14.12) 


w  is  the  complex  potential  of  ippendix  11.  T  is  evaluated  as  an 
integral  in  the  complex  z  plane  of  figure  17  and  with  (11.7)  of 
Ippendix  11  and  Gradsbteyn  and  Kyzhik  1980  (p307,  p942)  we  obtain 

|  -  i  (H/x) expl  L  (  iuH/xOc ) 1  si nh (oh/Oc )exp(-uH/Uc )  ,  h  <  H, 

i  i  (H/x) expl  L  (-iuH/xUc  )  J  sinh(oH/Uc  )exp<-<dh/U* )  ,  h  >  H  , 

( A4 . 13) 

In  the  limit  <a>H/Uc  >>  1.  when  the  duct  height  is  extremely  large 
relative  to  the  hydrodynamic  length  scale  Uc/o,  use  of  the  asymptotic 
form  of  L  (c.f  Gradsbteyn  and  Ryzhik  1980,  p940)  gives  for  u  >  0 


In  the  case  of  a  long  wavelength  sound  wave  radiating  iron  the  duct 
the  perturbation  potential  in  the  absence  of  displacement  waves  is 
and  in  the  vicinity  of  the  leading  edge  of  the  plate  it  can  be  treated 
as  the  potential  of  an  incompressible  flow.  The  complex  potential  of 
this  flow  will  be  of  the  form  pw(z),  where  w(z)  is  given  by  (Aid)  of 
Appendix  £1  and  represents  flow  into  the  duct  which  becomes  uniform 
with  speed  unity  well  within  the  duct.  p  is  chosen  to  match  the  near 
field  of  the  i,  velocity  of  the  incident  and  reflected  waves,  A.+A,  ol 
(3.9),  i.e  p  -  ^<A,*Ar)/ixi  in  the  limit  kt,x,  ^  0.  legleeting  order 
IP  relative  to  unity  we  find  from  (3.9),  (3.1)  of  the  main  text 


p  -  - 1  koAo  <1  -  BtJKJM  1 ) ) . 


(A4 . lb) 


(o>  is  found  from  condition  (A4.4)  and  (A4  7)  with  now  equal  to 
pw(z)  and  (c.f(A1.3)  of  Appendix  Al) 


dW*  |  -  -2pH/x  . 
dZ  Z=1 

We  find  for  J*( o> 

V(u>  =  (1  i)p(l  xl  H)~/2  , 

=  1(1  i)kofo(lxlH)“d  8**(HB))/2  (A4.16) 


7? 


